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We investigate the nonlinear self-trapping plienomenon of the Bose-Einstein condensates (BEC) 
in a symmetric double-well, emphasizing on its behind dynamical phase transition. With increasing 
the nonlinear parameter depicting the interaction between the degenerate atoms the BEC turns to 
be self-trapped manifesting an asymmetric distribution of the atomic density profile. Essence of 
this phenomenon is revealed to be a continuous phase transition and underlying critical behavior is 
studied analytically and found to follow a logarithm scaling-law. We then go beyond the mean field 
treatment and extend to discuss the effect of the many-body quantum fluctuation on the transition. 
It is found that the transition point is shifted and the scaling-law is broken. In particular, the 
quantum phase transition is accompanied by the change of the entanglement entropy which is found 
to reach maximum at transition point. Behind physics is revealed. 

PACS numbers: 03.75.Gg, 68.35.Rh 



I. INTRODUCTION 

Double-well system is a paradigm model used to 
demonstrate marvellous quantum tunnelling jj]- The re- 
alization of dilute Bose degenerate gas in last nineties 
provides a possibility of directly observing the tunnelling 
in the matter wave of macroscopic scale up to 100/im 
In Bose-Einstein condensates (BECs) system the inter- 
action between the degenerate ultra-cold atoms plays an 
crucial role. It dramatically affects the quantum dynam- 
ics and leads to many unusual phenomena like nonlinear 
Josephson oscillation, nonlinear quantum tunnelling and 
critical onset in coherent oscillations , etc. (3, 1^ ^ 1^- 
These problems have attracted much theoretical atten- 
tion over the past few years and the recent realization of 
the BECs in the optical trap of a double-well configura- 
tion has brought a new research surge 

Among many findings, the transition to self-trapping 
is most interesting one ITol Ull . It says, with in- 
creasing the atomic interaction (repulsive), the Josephson 
oscillation between two wells is blocked, the BECs atoms 
in a symmetric double-well potential shows highly asym- 
metric distribution as if most atoms are trapped in one 
well. This phenomenon was observed in lab recently @. 

In this paper we achieve insight into this somehow 
counterintuitive phenomenon with addressing its behind 
phase transition and the influence of many-body quan- 
tum fluctuation on the phase transition. Analytically 
we identify the self-trapping phenomenon as a continu- 
ous phase transition and critical behavior is found to be 
characterized by a logarithm scaling-law. We then extend 
to discuss many-body quantum fluctuation effect on self- 
trapping. We find that, the transition point is shifted and 
the scaling-law is broken down due to the quantum fluc- 
tuation. Further investigations show that quantum prop- 
erties are quite different for different interaction regions, 
such as self-trapping region and others. These properties 
can be well illustrated by entanglement entropy quan- 
titatively. With employing average entropy as the or- 



der parameter, we can clearly demonstrate the quantum 
phase transition: The entanglement entropy reaches its 
maximum at quantum transition point. 

Our paper is organized as follows. In Sec. II, the self- 
trapping is studied within mean-field framework both an- 
alytically and numerically, revealing the critical behavior 
at the transition point. In Sec. Ill, we discuss the many- 
body quantum fluctuation effect and reveal the quan- 
tum entanglement manifestation of transition to nonlin- 
ear self-trapping. Sec. IV is our discussion and conclusion. 

II. TRANSITION TO SELF-TRAPPING AND 
SCALING LAW 

For two weakly coupled BECs trapped in a symmetric 
double-well, the system can be described by the so-called 
two- mode Hamiltonian [sLIl^l 

H = l (st a - Vb) + ^ (st « - b^b) ' I (st 6 + b^a) , 

(1) 

where the Bose operators a^^^ and 6^^^^ correspond to 
annihilating (creating) operators for different well re- 
spectively, 7 = E^^ ~ E° is the energy bias between 

the two wells and Ef = / (l^lVi^ip + V {r)\(pi\'^'^ dr , 

c — Ci— J |(^i|^dr denotes the effective interaction 

of atoms, V — J {^^^'Pa^'Pb + V{r)ipaVbj dr is the ef- 
fective Rabi frequency which describes the coupling be- 
tween two wells, N is total atoms number which is con- 
served, a is the s-wave scattering length, and ipi {i = a, b) 
are wave functions for each well respectively. In the 
present work we focus on the case which has been realized 
in lab recently. For this case, the potential is symmetric 
so that 7 = 0, and the interaction is repulsive, i.e., c > 0. 

If the particle number is larger enough, the system 
can be well described in the mean-field approximation. 
Under mean-field approximation, the dynamics of the 
system is described by a classical Hamiltonian H = 



2 



H /N (up to a trivial constant) in which 
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GP/ 
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|0) is collective state of A^- 

particle system ^ = |a|e'^" and b — |5| e**"" are 

two c-numbers which correspond to the probability am- 
plitudes of atoms in two wells respectively. By intro- 
ducing the population difference s = \b\^ — |o|^ and the 
relative phase 9 = db — da, the classical Hamiltonian can 
be reduced to 



v\/ 1 — cos 6*, 



(2) 



where s and are canonical conjugate coordinates. Their 
equations of motions are 



vs 



-cs - 



■,e. (3) 



Self-trapping motion refers to the trajectories whose 
average population difference is not zero < s 0. In 
the experiment jS], all the atoms are placed initially in 
one well, i.e., s(0) = I or —1. As observed, with small in- 
teraction, the Josephson oscillation will be observed, and 
with larger interaction, the self-trapping emerges. This 
phenomenon can be well understood by the above classi- 
cal Hamiltonian systems Fig.l plots the evolution of 
population difference s and its average for different inter- 
actions calculated by lj2Jl with initial condition s(0) = 1. 
For c/v is smaller than 2, the population difference is os- 
cillating symmetrically between 1 and —1 and its average 
is zero. However, for c/v is larger than 2, the motion is 
limited in half plane and the amplitude decreases with the 
interaction increasing, hence, the average of population 
difference will be nonzero and increasing with interaction. 
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FIG. 1: (Color online) For initial condition s — 1, the popula- 
tion difference evolves with time for c/v = 1, 1.995, 2.005, 2.5 
(left column), and the average of population difference vies 
parameter c/v (right column). Inset: the critical behavior 
near c/v — 2, the red line is for analytic formula (5) and the 
circles are for numerical simulation 

The above process can be well understood from the 
analysis on the phase space of the classical Hamilto- 
nian system. In Fig. 2, we plot the trajectories in phase 
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FIG. 2: (Color online) Trajectories on the phase space of the 
classical Hamiltonian system (2) (upper panels). In bottom 
panels we plot the energy profiles for the relative phase 6 — 
(dashed) and 6 = it (solid), respectively. The energies of the 
trajectories in upper panels are also denoted in bottom panels 
respectively using the lines with the same colors. 



space and classical energy profiles for different parame- 
ters. The red lines correspond to the trajectories of which 
all the atom are initially in one well, i.e., s(0) = 1 or 
— 1. From this figure, we can see clearly that the dynam- 
ics transition happens at the moment when the energy 
of the trajectory with initial condition s(0) = 1 or — 1, 
H{s — ±1) = c/2, equals to the energy of the separa- 
trix which \s H = —v. When the energy of the trajectory 
initially in one well is larger than —v, one only finds 
Josephson oscillation trajectory, when for the energy be- 
ing smaller then —v, the self-trapping happens. 

For classical Hamiltonian system we can obtain the pe- 



riod T of a given trajectory by the integral T — f -gjjds, 
and average s of it by (s) — f -gjjsds, in which the inte- 
gral path is along the trajectory. For the trajectory with 
initial condition s(0) = ±1, we have H{s = ±1) = c/2. 
Thus, from 101 and Q, we get 



T = 



and 



ds 



!)^(l-s2)-[c(l-s2)/2t;]2 
9 ds 



-s2)-[c(l-s2)/2t,]2 



c/v < 2 
c/v > 2 



(4) 



2_ fl 

T J-1 



c/v < 2 



Tl^l-{2v/cy „^(i_,2)_y(i_,2)/2„]2 C/V>2 

(5) 

in which we have used the formula cos9 — c\J\ — s'^ /2v. 
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After some elaboration, we obtain 

[ cjv < 2 

[ 2c/vIrn{K[(c/v)'y((c/vy^-i)]) > ^ 

where K{x) is the complete first kind of elliptic integral. 
Near the transition point, it exhibits the logarithmic crit- 
ical behavior 



± 



V27r 



c/i;ln[(c/w)2 - 4] ' 



(7) 



The inset figure of Fig. 1 plots this critical behavior, 
where theoretical result is confirmed by numerical result 
obtained by numerically solving Eq.(3) with 4th-5th step- 
adaptive Runge-Kutta algorithm. 

Our logarithmic critical behavior is very similar to the 
critical behavior in the measure synchronization in cou- 
pled Hamiltonian systems 0. This is because critical 
behavior in both case is closely related to the separatrix 
of the Hamiltonian. Near the separatrix the period of the 
trajectory diverges to infinity, as the function of the rel- 
ative deviation of the energy from the separatrix energy, 
its divergency follows a logarithm lawjl^. 

In the above discussion, the initial state is set as s = 1, 
in fact for any initial state denoted by Si , 9i the transition 
to self-trapping occurs at some interaction parameter, 
and the critical behavior follows the same logarithm. If 
we extend the above discussion to this general case, we 
can obtain the general criterion for the occurrence of the 
self-trapping, i.e., H{si,9i,c,v) < —v. Then, the critical 
point is expressed as. 



(-). 

V 



2(1 



1 



(8) 



From the above analytic expression, we see that, for the 
initial state with smaller population difference it requires 
stronger nonlinearity so that self-trapping occurs. More- 
over, the critical point can be adjusted by the relative 
phase between the two weakly linked BEG in double- 
well. For example, for the case the population difference 
is 0.5, the critical point approximates to 15 and 8 for 
Qi — Q and 7r/2, respectively. In practical experiments, 
the relative phase can be adjusted with using 'phase- 
imprinting' method, i.e., shedding un- uniform laser light 
on the BECs in double-well. This method has been suc- 
cessively applied to generate the dark solitons in cigar- 
shape BECS(17j. 



III. MANY BODY QUANTUM FLUCTUATION 

EFFECTS 



In the mean field treatment we assume that the number 
of particle is large enough. However, in practical experi- 
ment, the particle number is finite, in order to know the 
quantum fluctuation effect due to finite particle number, 
we should investigate the self-trapping within the frame- 
work of the many body quantum system [Tsl IT^ . 



In treating the quantum many-body problem it is 
helpful to bear in mind some results from the quan- 
tum information theory concerning the entanglement 
|20|. The quantum entanglement is realized not only 
to be a crucial resource that allows for powerful com- 
munication and computational tasks that are not pos- 
sible classically, but also to be a signal for quantum 
long-range correlation and therefore can serve as in- 
dicator for the quantum transition in concrete solid 
system 21]. Recent years witness growing interests in 
studying the interplay between entanglement and quan- 
tum phase transitionlil I2ll2ll2ll2q. 

Previously, some eff'orts have been devoted to study 
the dynamics of the BECs in double-well with full quan- 
tum treatment [271 12^ . In Ref . [23| the authors presented 
a quantum phase-space model of the BECs in a double- 
well potential by using the Husimi distribution function. 
They showed a good correspondence between the phase 
space of classical Hamiltonian (mean-field approxima- 
tion) and the quantum phase-space of two-mode hamilto- 
nian J^l (full quantum framework). The authors of Ref. 
Us] calculated the time evolutions of states and their cor- 
responding entanglement with different initial states for 
several different interactions between atoms. The time 
evolutions of entanglement entropy presented in Ref. 
for several interactions between atoms show a decreasing 
tendency with increasing interactions. 

In our following discussions of this section, however, 
we focus on the critical behavior at the transition to self- 
trapping as revealed by the above discussions, addressing 
how the quantum fluctuation influence on the transition 
behavior. As will be shown latter, the transition point 
is shifted and the scaling-law is broken down due to the 
quantum fluctuation. With increasing the atom number, 
the transition behavior demonstrates a perfect classical 
quantum correspondence. We also calculating the en- 
tanglement entropy achieving insight into the quantum 
transition. Our calculations on the time evolutions of en- 
tanglement entropy conflrm the results of |28l |. and our 
further calculations on the time averaged entanglement 
strongly suggest that the entanglement entropy of this 
system serve as a good order parameter to describe such 
quantum transition. 



A. Quantum Phase Transition 

In the quantum framework the evolution of the system 
is governed by the Schrodinger equation 



(9) 



where |V'(0) 

(«^) [b') |0) (n = 0,... ,7V) are Fock 
states, and a„ are the probability amplitudes respec- 
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tively. Hence, the population difference is given by 



E 



|a„r(jV~2n) 
N 



(10) 



We choose |0,iV) as the initial state in the full quan- 
tum framework, which corresponds to s = 1 in the mean 
field model. In Fig. 3, we plot the average population 
difference calculated from the above Schrodinger equa- 
tion for different total particle numbers. From the above 
calculation, we find that, the quantum fluctuation has 
two significant effects on the transition to self-trapping. 
First, the critical point is shifted to the left-hand side due 
to the finiteness of the particle number. From Fig. 3, it 
is clearly seen that, keeping parameters c,v as constant, 
with increasing the number of the atoms, the transition 
point clearly shift to left-hand. For example, for N — 50, 
the transition point shifts to c/v = 1.6. We know that 
the quantum fluctuation closely relates to effective planck 
constant, here for this model it is Therefore we 

expect the deviation from the mean field critical point 
should be inversely proportional to total atom number. 
This prediction is confirmed by our calculations as shown 
in the inset of Fig. 3. in which i5 is the difference between 
the quantum transition point and the mean field one. 

Secondly, the logarithm scaling-law is broken down by 
the quantum fluctuation. In Fig. 3 we flnd the quan- 
tum fluctuation destroys the logarithm scaling-law of the 
mean field and no clear scaling-law is observed for the 
quantum case. With increasing the atom number the 
quantum results tends to the mean field results in the 
limit iV to cxD, as required by the classical quantum cor- 
respondence principle p^. 

We should address that in the above calculation time 
period for average should be much longer than the period 
of the fastest oscillations but shorter than the period of 
the shortest quantum beating. This is because, essen- 
tially, the dynamics of quantum system is periodic or 
quasi-periodic, therefore, any dynamical effects of quan- 
tum system depends on the time scales 0, HO] • In our 
problem, there are two time scale, one is for integrating 
classical equation (3), the other is for integrating quan- 
tum equation (10), the later is N times the former one. 
In our calculations, the average time is 50A^, meaning 
that the corresponding classical time scale is same (=50). 
On the other hand, we find that, the averaged population 
difference increases gradually to 0.001 and then soars up. 
This is different from the mean field situation, where the 
averaged population difference keeps zero and then turns 
to be nonzero after critical point. So, in the quantum 
case, we define the transition point as the point that the 
averaged population difference is larger than 0.001. This 
observation also suggests that the no scaling-law for the 
finite particle situation. 



A 
(/) 
V 



1.0 



0.8- 



0.6 



0.4- 



0.2 



0.0 



■- ▼-- 


N=50 




N=100 


. A- ■ 


N=200 


■- ■-- 


N=400 




Mean-field 




50 100 150 20D 250 3110 350 400 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



c/v 



FIG. 3: (Color online) The average of population difference 
vies the parameter c/v obtained by full quantum simulations. 
Here, the initial state is (n = 0, • ■ ■ ,N), namely s = 1. In- 
sert: the shift of the critical point comparing with meanfield 
predication S for different atom number. 



B. Entanglement Manifestation of the Quantum 
Phase Transition 

To well understand the self-trapping phenomenon in 
full quantum description, we calculate evolution of the 
occupations on Fock states \n,N — n) [n = 0, • • • ,N). 
Fig. 4 shows the evolution of Fock state occupation for 
different interactions. The horizontal axis is the time i, 
the vertical axis is the index of Fock state, namely, n cor- 
responds to |n, iV — n), and the contour is for the occu- 
pation probability. Fig. 4(a) is for the linear case, which 
shows the occupations are oscillating between |0,iV) to 
|iV, 0). For c = 1.95 (seeing Fig. 4(b), near the transi- 
tion point, we see that the wavefunction spread much 
rapidly to all the Fock states. Fig. 4(c) is plotted for the 
self-trapping case, c = 2.5, from which we see that the 
occupations is narrowed in partial Fock states, so that 
the average of population difference is nonzero. 

From Fig. 4, we also see that the dynamics properties 
of such a quantum system are quite different for differ- 
ent interaction regions. To achieve more insight into the 
quantum transition to the self-trapping, we introduce the 
quantum entanglement entropy. For the system with the 
wave function \ip) = X]^=o 1'^' ~ density op- 

erator is given by 



(i^l =^a„a;,|n,iV-7i) (7V-m,m|, (11) 



Taking the partial trace with respect to one well yields 
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FIG. 4: (Color online) Occupations on Fock states for differ- 
ent time. The contours denote the probabilities of occupa- 
tions. From up to bottom, cjv = 0, 1.995, 2.5 respectively. 

the reduced density operator for the other, 

Pa-5]|a„|'|n)(n|. (12) 

n 

Thus, the entropy of entanglement between the two cou- 
pled EEC's is given by il9|| 

TV 
n=0 

The entanglement entropy has the following properties: 
its reaches maximum Ei^p) = log A'^ when |a„p — and 
its minimum E{p) ~ when |a„p = 1 and others are 
zero. 

Because the self-trapping is dynamic phenomenon, the 
occupation on each Fock state is varied in time, there- 
fore, we use average entropy. Technically, we have two 
choices to average it: we can average occupation firstly 
and then calculate entropy, and calculate entropy firstly 
then average it. So, we denote 

AT 

= -^(|a„P)log(|a„p)/logiV, (14) 

and 

{E) = /-^|a„plog|a„|A/log7V. (15) 

\ n=0 / 

The above formula have been normalized by log N . 



FIG. 5: (Color online) The entanglement entropy for different 
cjv. The circle line is for Ea-u and dashed line for {E) 

In Fig. 5, we plot the two kinds of average entropy 
with different interaction. Obviously, the dynamics of 
the quantum system can be well illustrated by the av- 
erage entropies in quantities. For the linear case, the 
atoms mainly occupy several Fock states at a given time 
(see Fig. 4 (a)), and the occupied states are changed with 
time. Thus, the instantaneous entropy is small, so does 
the average entropy {E) . On the other hand, because the 
average population on each Fock state are almost equiva- 
lent for this case, so E^v should be large. With the inter- 
action increasing, the occupations on Fock states extend 
so that the instantaneous entropy increases. However, 
when the interaction exceeds the transition point self- 
trapping occurs and the occupations are limited on sev- 
eral Fock states (see Fig. 4 (c)), hence, the instantaneous 
entropy becomes small, so does (£') . For the same reason 
Eav will be small for self-trapping cases. From Fig. 5, 
we also observe that Eav is almost independence of the 
particle number and varies very little with changing the 
interaction parameter before it reaches its maximum. Af- 
ter that, it shows quite sensitive on the particle number 
as well as the interaction parameter. Whereas, for the 
entropy < E > before and after the maximum point it 
shows strong dependence of the particle number. 

It is interesting that the two average entropies reach 
their maximum at the point very close to the transition 
point of mean field c/v = 2. This is very similar to 
phase transition of spin systems, where the phase transi- 
tion happens at thepoint when entanglement of system 
reaches maximum '24l|. This is different to the situation 
of the two-impurity Kando model where the entangle- 
ment vanishes at a quantum critical point [25^. In par- 
ticular, we find that, the maximum points of the aver- 
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age entropies varies very little with changing the particle 
number. The property suggests that the entanglement 
entropy is a better quantity than the average population 
in serving as a indicator to quantum phase transition, 
because the latter is too sensitive on the particle number 
as shown in Fig. 3. 

IV. DISCUSSION AND CONCLUSION 

In conclusion, we have made thoroughly analysis on the 
transition to self-trapping for BECs confined in a sym- 
metric double-well. Analytically we identify it as a con- 
tinuous phase transition, where the time averaged popu- 
lation difference between two wells changes from zero to 
nonzero following a logarithm law at a critical point. We 
also discuss influence of the many-body quantum fluc- 
tuation on the transition to self-trapping. We find that 
the transition point is significantly shifted by quantum 
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